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Abstract 

The sawtooth chain with pairs of S = 1/2 spins interacting with XXZ-interactions placed on each second tooth 
is considered. All other interaction bonds are taken to be of Ising type. Exact statistical mechanical solution of 
the model within the direct transfer-matrix technique is obtained. The solution allows one to obtain exact analytic 
expressions for all thermodynamic functions of the model. Ground stated properties are also investigated, the 
corresponding ground state phase diagram is presented. 

"TT 1 : 1 Introduction. 

The sawtooth chain or delta chain is a one-dimensional lattice spin system with a topology of corner-sharing triangles 
(Fig. [I]). This system is famous for a number of important features. Physically, magnetic lattices corresponding to 
the sawtooth chain are found in a number of compounds, in delafossite YCUO2.5 [TJ H] and olivines with structures 
Z11L2S4 (L=Er,Tm,Yb) [3] to cite a couple of examples. Antiferromagnetic Heisenberg model on sawtooth chain is 
strongly frustrated, however, in contrast to many other frustrated systems, the ground state of sawtooth chain is 
known exactly. Namely, the ground state in case of cyclic boundary conditions is a two-fold degenerated dimcrised 
structure where either each left pair of spin at each triangle or each right pair forms spin-singlet states (dimers) [1] . The 
excitations of the sawtooth chain are gapped and have topological origin. They are found to be "kink" -"antikink"- 
type domain wall structures [5j [6]. Anisotropic variants of sawtooth chain have also been investigated revealing 
additional features [3 [8j . Another important feature of the sawtooth chain as well as some other frustrated spin and 
Hubbard lattices revealed quite recently is the appearance of localized magnons states (dispersionless excitation bands) 
[HI EH H3 [H] ■ The localized low energy excitations in sawtooth chain which exist due to frustrated geometry, more 
precisely, due to the triangular plaquette affect the low-temperature thermodynamics. The corresponding eigenstates 
have been constructed for sawtooth chain in case of Heisenberg model in Ref. [T2] and for Hubbard electrons in Ref. 
[TP] , Calculating the ground state degeneracies for the localized magnon (or electron) states one can obtain explicit 
expressions for thermodynamic quantities for low temperatures and near saturation field [101112] . However, the problem 
of exact description of whole thermodynamics for sawtooth chain as well as for many other interacting quantum spin 
systems is still an open issue. 

One can mention a formal approach to these problem which is not justified properly yet but in some particular 
cases demonstrates rather good agreement with experimental data and numerical calculations. The approximation 
consists in replacement of some or all interaction bonds with Ising ones |13j-|24j. As a result one can obtain an 
interacting spin system which allows one to calculate all thermodynamic functions analytically. In this paper we use 
this approach for exact solution of the system with Ising and Heisenberg interaction bonds on sawtooth chain. To make 
possible exact calculation of the partition function within transfer-matrix method we consider all interaction bonds of 
sawtooth chain as Ising ones except left bonds on each second triangle. Spins connected by these bonds interact by 
Heisenberg XX Z interaction. Such a construction allows one to represent the Hamiltonian in the form of sum of block 
Hamiltonian commutating to each other, thus, the exponent in partition function can be expanded to the product 
of exponents corresponding to each block. Due to this fact, the method of transfer-matrix can be exploited leading 
to exact calculation of all thermodynamic functions. In the first section we formulate the model and give its exact 
solution, then, in the next section, we analyze various ground states of the system and draw a corresponding phase 
diagram. The last section contains concluding remarks. 



2 The model and its exact solution. 

Let us consider the sawtooth chain in which the left pair of spins on eevery second triangle interacts with the XX Z - 
Heisenberg interaction, while all other interaction bonds are of Ising type, i.e. the interactions include only z-component 
of spin operators. The corresponding Hamiltonian is suitable to write as a sum of block Hamiltonians where each block 
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contains 5 spins: a pair of spins connected with quantum bond, two adjacent spins from the basement and one spin at 
the top of the triangle from the left hand side (See Fig. \2§ : 

N / H \ 
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Hi = J (A (5f x 5f 2 + S^5&) + SM) + K ^ fa + <r i+ i) + l^r, + if 4 S>i +1 - + Sf 2 ), 



where S^, a = 1, 2, a = x, y, z are S = 1/2 spin operator components of the spin connected by quantum bond from 
i — th block, Ui and Oi+\ are two adjacent spins from the basement and r, is the spin on the top of left triangle; 
as only z-components of all a and r spins are included in the interaction, one can consider them just as a classical 
variable, taking values ±1/2. Here we also choose the most general form of interaction with different coupling constants 
J, K\, K2, -K3, K4, for different kind of bonds and different g-factors for Zeeman term corresponding to different kind of 
spins. The partition function of the system then reads: 

^ = EE s ps^- ( 2 ) 

W (t) 

Here one should sum over all values of classical variables a and r and take a trace over all states of spin operators S. 
Because of commutativity of Hamiltonians for different blocks one can expand the exponent and get the product of 
the terms corresponding to different blocks. After that the traces for each blocks can be taken separately: 

N 

Z = EII E ^(^^ + i\n)e^-^' + ^' + ^\ (3) 

(a) i=l t*=±1/2 

where 

Q(a i ,a i+1 \r i ) = Spe-^ = J2e- l3K{cri '' Ti+lln \ (4) 

71=1 

and A„ are four eigenvalues of the TLi. 

Ai, 2 fa, <Ti+i\n) = 2 J =F ~ (2i? - Kx fa + <j i+ i) - K 3 Ti - K 4 a i+l ) , (5) 
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A 3 , 4 fa,0"i+l \n) = —-J =F 2 V (^1 fa + °"i+l) + K STi - K4<T l+1 ) + J 2 A 2 . 

In Eq. ([3]) the sum over all states of spin r in each block can be taken independently from the other, then yielding 

N 

2=^n^( ff ..%^^ ( ' i+ ' i+i) . (6) 
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Z fa, <r m ) = n fa, <r m |l/2) e ^ H ^ K ^ + ft fa, - 1/2) e -P^-^). 

Thus, the partition function takes the form similar to that of the single chain with N sites with classical variables <Ji 
in them. The partition fuction can be calculated within the standard transfer-matrix technique (See for example 25j): 

N 

z = e n t . ^ ) = TrTAr = a + + a - (7) 

(a) i=l 

here the cyclic boundary conditions fJi+w = are imposed and 2 by 2 transfer-matrix T with eigenvalues A± has the 
following form: 



Z(l/2,l/2)e^ 2(1/2,-1/2) 
£(-1/2,1/2) Z(-l/2,-l/2)e-^ 



(8) 



Thus, for free energy per one block of the sawtooth lattice with mixed Heisenberg and Ising bonds described above in 
the thermodynamic limit when only maximal eigenvalue survives from Eq. ([7]) one obtains 



/ = — log l 1 - lz + eP* + Z-e~& +J(z + eP* Z^^)' + AZ Z Q ) , (9) 



2 



where the following notation arc introduced: Z + = £(1/2,1/2), Z- = Z(-l/2, -1/2), Z = Z(l/2,-l/2) and 
Zo = Z(— 1/2, 1/2). Then, one can obtain analytic expressions for all thermodynamic quantities of the systems by 
taking derivatives of Eq. ©. Thus, Eqs. flS])-© exactly solve the problem of thermodynamics of the systems under 
consideration. Let us finally list the explicit expressions for the entries of transfer-matrix for particular choice of 
couplings, mentioned in the beginning of the next section, Ki = J\ and K 2 = K 3 = K4 = J and for the same value of 
g- factor for all spin, which imply Hi = H 2 - 

Z+=2 f e H(H-i/2J) cosh(/? ^ H _ 1/2(Ji _ 1/2 J))) + P H cosh \8l/2\J J% + J 2 A 2 \\ (10) 

+e -Pi(H-i/2J) f e -pi CQsh ( H _ 1/2(Ji _ 3/2J))) + e p i cosh (/31/2V( Ji - J) 2 + J 2 A 2 ))) , 
Z_ = 2 L?k&+V**) (e~ i cosh(/3 (H + l/2(Ji - 1/2J))) + e?* cosh (/31/2v/(./i - </) 2 + J 2 A 2 
+e -/3i(H+i/ 2 j) / cosh(/? ^ + 1/2(Ji + 3/2 j))) + e /3i C osh \8l/2^Jl + J*&yyj , 
Z = 2 ( e /3|(ff-V2./) r e -/?i cosh(/3 _ 1/4 J))) + e p i cosh (j31/2J^/l + A 2 ^ 
+e -/J§(ff-i/2J) ^ e -/Ji cosh ^ + z/AJ)) + e i cosh (01/2 JA))) , 

Z = 2 ( e /3 ^ (ff+1/2 ' 7) (e"' 3 * cosh(/3(# - 3/4 J))) + cosh (/31/2 JA)^ 
+e -/3i(if+i/ 2 j) /-^ CQsh ^ (ij + 1/4J)) +e /3f cosh (/31/2JVTT 

3 Ground states phase diagram for antiferromagnetic couplings 

Let us analyze different ground states of the model. For the sake of simplicity we restrict ourselves to the case of only 
two different coupling constants, one along the basement J\ and another one between spins on the top and spins on the 
basement J. This means one should put K\ — J\ and K 2 = K3 = K4 = J. We also consider only antiferromagnetic 
coupling J > 0, J 1 > 0, the rest cases are to some extent trivial ones. Solving the eigenvalues and eigenvectors problem 
for block Hamiltonian Tii and taking into account all combinations of spins cr^ , cr t ; + i and Ti one can get the following 
periodic eigenstates for whole chain with period equal to the period of the lattice up to the inversion of all spins: three 
antiferromagnetic states with M — 0, two of which are degenerated and differ from each other by the flip of all a and 
r spins (the corresponding energies per one block are presented) 



\AM+) = [] \S Z = 0, +)i (g) \cn =t, n =[) (11) 

»=i 

N 

|AM_) H = \S Z = 0, -) 4 (g) [at =i,n =t), 
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CAM = ~\ \J + ^{J-Jxf + J 2 ^ 

where \S Z = 0, ±)j stands for the following state of two S spins form z-th block with total projection equal to zero 

\S Z = 0, ±) = j± (I TI) - 7 ±l IT)) , (12) 

2 T9A2 



±(j-j 1 ) + j(j-j 1 y + .pA 

7± = JA ' 

This eigenstate corresponds to the A3 eigenvalue from Eq. ^ . For arbitrary values of a and r surrounding pairs of S 
spins coefficient 7 from Eq. (fT2"|) is 



J (o-i+i - n) - Ji (pi + a i+ i) + J( J (a i+1 - n) - J l ((Ji + cr i+ i)) 2 + J 2 A 2 

l{<Ti,<7i+l\Ti) = J£ , (13) 

Thus, in these configurations the pair of spins connected with the Ising bonds are aligned opposite to each other, 
while each pair of spins connected with the Heisenberg bond is in the special state with S* ot = 0, see Fig. ([3]). 
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Another antiferromagnetic ground state consists of triplet spin configurations on each Heisenberg bond and pairs of 
spins connected with Ising bond both pointed opposite to the direction of S = 1 pairs. However, this ground state 
never realize at T = for antiferromagnetic region of values of J and J\ . There are also three ferrimagnetic state of 
the chain when external magnetic field is turned on. All these states have spatial period corresponding to the period 
of the chain and magnetization equal to 1/2. In the first ferrimagnetic state all pairs of S spins are in a singlet state, 
while the rest a and r spins are pointed along the field(Fig. 2]): 

N 1 / 

|F1) = J] |S Z = 0), ®k(=t,Ti=T>, e F1 = --^/j! + J 2 A 2 -ff, (14) 
»=i 

where eigenstate \S Z = 0)i is given by the same formula as Eq. (fT2"|) but with another value of the coefficient gamma, 
which is given by Eq. (TiU|) 



7= z*±vg±Z_ . (15) 

In the rest two ferrimagnetic states all pairs of S spins are in the triplet polarized state with S z = 1, while a and 
t spins surrounding them are pointed opposite to each other: 



|i f, 2)=niTT>i®ki=t,Ti=J.), e F 2 = \ji-H (16) 

i=l 

N 1 

1^3) = [] I TT>« <S> =^ T * =T). = —yh - H. 



And finally in the spin polarized phase all spins are pointed along the field, 

N 



|5P) = JJ|Tt><(g)k<=T,7i=T), esp = J+\jx-2H. (17) 



i=l 



Besides the states with non-broken lattice symmetry the system under consideration can posses another ground state 
with a spatial period equal to two blocks. This state can be realized at T = only for i] = ^j- < 1. This state 
can be obtained from \AM_) by flipping every second a spin in the basement of the chain (Fig. [5]). This state with 
broken translational symmetry corresponds to magnetization of the system equal to 1/4. We will refer to it as a spin 
modulated phase, \SM). The spacial period of the state and value of magnetization per site are in full agreement with 
Oshikawa-Yamanaka- Affleck criterion [26] . The corresponding energy per one block is 

-ijfl + IAI + v/I+A^-itf. (18) 



esu 4 - v , 2 

Now it is straightforward to draw the T = ground state phase diagram for antiferromagnetic couplings J > 0, J\ > 0. 
In Fig. © one can see the phase boundaries in (rj,h) plane (h = H/J) for A = 0.3. An interesting feature of the 
system is that at antiferromagnetic values of couplings the ground state for H = is unaffected by the value of rj. 
Thus, the unique ground state is either |AM_) or |AM+). One should distinguish two regions of r], namely, for r\ < 1 
increasing the magnetic field one can see transition to spin modulated phase with M =1/4 prior to appearance of the 
ferrimagnetic structure Fl with M = 1/2. However, for rj > 1 antifcrromagnctic phase changes immediately with the 
ferrimagnetic one, thought width of the AM phase region (the value of the h at which the zero temperature quantum 
phase transition takes the place) with rj — > oo asymptotically goes to zero. So, one can see a small closed region of 
SM phase in the (rj,h) plane. The values 77 = 1, h = the system demonstrate additional feature, the macroscopic 
degeneracy of the ground state. At this point the the S spins in each block are in the singlet state, whereas its adjacent 
a and r can take arbitrary value with one restriction a + r = 0. Thus, in each block a and r can freely take either 
a = 1/2, r = —1/2 or a = — 1/2, t = 1/2 values. Thus, one obtains macroscopic two-fold degeneracy at J = J%. One 
can refer to this state as a frustrated antiferromagnetic state: 

N 

\FR) = H\S*=0) i (g)\S i }, (19) 
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where can be either \a- L =T,Ti =J,) or \a- L =l,Ti =|) And finally the equations of phase boundaries between four 
phases presented in the phase diagram in Fig. §5§ are listed below. 



between \AM) and \SM) 



between \AM) and \Fl) : 



between \SM) and \F1) : 



between \Fl) and \SP) : 

4 Conclusion 

In this paper using the example of sawtooth chain we considered the possibility to obtain exact thermodynamic solution 
for the one-dimensional quantum spin systems by considering their counterparts where some interaction bonds are 
changed with the Ising ones. Namely, we considered the pairs on S = 1/2 quantum spins interacting with XX Z 
interaction arranged in the sawtooth chain in such a way that allows one to use transfer-matrix technique for exact 
calculations. This result continues the series of investigations of the subject performed earlier for other one-dimensional 
spin systems with Ising and Heisenberg bonds |1.3j- [24] . These results are not only of academic interest. There are many 
one-dimensional spin systems which are the models of real magnetic material and which are not integrable. Until now, 
only the actually reliable way to describe the thermodynamic and magnetic properties of such a system is laborious 
numerical calculations. However, results of recent investigations are evidence that changing some bonds of the system 
with Ising ones, on the one hand, makes the thermodynamical problem exactly solvable and, on the other hand, does 
not change drastically magnetic properties of the system at least in case when only ferromagnetic bounds have been 
changed. One can look at the plots of low-temperature magnetization curves for the system considered in the paper 
(Figs. 0, ©). As it can be seen from phase diagram (Fig. |SJ), there is qualitative difference between rj < and 
r\ > 1 cases. For r\ < 1 one obtains the curve with three intermediate magnetization plateaus, at M = 0, M = 1/4 and 
M = 1/2, which corresponds to the stability regions of AM, SM and Fl structures respectively. In Fig. one can 
see the corresponding plot for A = 0.3 and r\ = 0.5. For r\ = 1 and A = 0.3 one can see in Fig. ((8} only one plateau at 
M = 1/2. This result qualitatively corresponds to that obtained in purely quantum sawtooth chain [12] . Numerically 
obtained magnetization curves for sawtooth Heisenberg chain contain only plateaus at M = and M = 1/2 for whole 
range of parameters. Though, there is rather large discrepancy in quantitative characteristics, such as positions of the 
terminal points of the plateau, but this can be explained by the fact that in our case we changed antiferromagnetic 
bonds with Ising ones, which generally speaking makes much more changes into the thermodynamics of the system 
than changing ferromagnetic bonds. An additional feature which is absent in case of purely quantum antiferromagnetic 
sawtooth chain is the plateau at M — 1/4 for r) < 1. In order to obtain more satisfactory correspondence between 
numerical results for purely quantum system and its mixed counterpart, one can consider a larger quantum cluster, 
containing three or more spins. 
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Figure 1: The sawtooth chain. 




Figure 2: The sawtooth chain with Heiscnbcrg bond on every second triangle. The Heisenberg bond is marked with 
bold line. Group of spins included into the dotted rectangle correspond to one block with Hamiltonian Hi 




Figure 3: Spin configuration corresponding to the \AMJ) antiferromagnetic state. Ovals mark the pairs of spin 
connected with Heisenberg bonds in \S Z = 0, — } state. One should inverse all arrows to obtain |5 Z = 0, +) state. 




Figure 4: Spin configuration corresponding to the \F1) ferrimagnetic state. 




Figure 5: Spin configuration for corresponding to \SM) phase 
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Figure 6: T = ground state phase diagram for the sawtooth chain with Ising and Heisenberg bonds for J > 0, J\ > 
in (i], /i)-plane for A = 0.3. 
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Figure 7: Magnetization process for rj = 0.5 at extremely low temperature T = 0.001J and A = 0.3. Three magneti- 
zation plateaus at M = 0, M = 1/4 and M = 1/2 correspond to AM, SM and Fl phases. 
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Figure 8: Magnetization process for 77 = 1 at extremely low temperature T — 0.001J and A = 0.3. Spin modulated 
phase corresponding to the plateau at M — 1/4 is absents here. One can see only two plateaus at M = and M = 1/2 
corresponding to FR and Fl phases. 



9 



